Electronic structure and non-magnetic character of 5-Pu-Am alloys 
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The around-mean-field LSDA+U correlated band theory is applied to investigate the electronic 
and magnetic structure of /cc-Pu-Am alloys. Despite a lattice expansion caused by the Am atoms, 
neither tendency to 5/ localization nor formation of local magnetic moments on Pu atoms in Pu- 
Am alloys are found. The 5/-manifolds in the alloys are calculated being very similar to a simple 
weighted superposition of elemental Pu and Am 5/-states. 
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Electron-electron correlation effects in a proximity to 
the localization threshold of the 5/ series, which is 
crossed between Pu and Am, recently attracted signif- 
icant attention in condensed matter physics 0, U • 
While strongly correlated systems usually have local 
magnetic moments, the whole set of experimental data 
clearly demonstrates their absence in Pu Q. Recent 
LSDA+U theoretical calculation have explained the 
non-magnetic character of <5-Pu, leading to the 5/ 6 -like 
state, hybridized with conduction-electron states. An im- 
portant issue remains open whether a further 5/ local- 
ization or often speculated onset of magnetism can take 
place when expanding the lattice by Am doping. 

Alloying Pu with Am provides an important opportu- 
nity to stabilize the high-temperature (5-phasc down to 
T = 0. Since Am atoms are larger than Pu, they create 
a "negative" pressure increasing the Pu interatomic spac- 
ing. Naively, knowing the sensitivity of the Pu phases to 
the lattice spacing, one can expect the formation of local 
magnetic moments at Pu due to an increase of 5 /-states 
localization. In contrast, recent susceptibility measure- 
ments of the Pu-Am systems do not show any systematic 
variations and susceptibility remains very weakly tem- 
perature dependent, not cxceedingLKU 8 m 3 /mol 0. 
Also results of specific-heat studies |g indicate no sign of 
magnetism in Pu-Am alloys. 

The Pu-Am alloys were studied theoretically in Ref. 
using the local spin density (LSDA) and generalized gra- 
dient (GGA) approximations to the Density Functional 
theory (DFT), and no electron-electron correlations be- 
yond those incorporated in LSDA/GGA were included. 
Large local magnetic moments on Pu and Am atoms in 
Pu-Am alloys (disordered or antiferromagnetically cou- 
pled) were calculated in Ref. , in contradiction to ex- 
perimental findings ■ The situation is similar to the 
DFT calculations (both within LSDA and GGA approxi- 
mations) for elemental <5-Pu 0,0] and fee- Am fl3| . 
where magnetism tends to appear as an undesirable ar- 
tifact. 

In this Communication we go beyond the conventional 



LSDA/GGA approximations and apply the LSDA+U 
correlated energy-band approach to investigate the elec- 
tronic structure of the Pu-Am alloys. The LSDA+U 
method 0,0] starts from the LSDA total energy, which 
is supplemented by an additional intra-atomic Coulomb 
U correlation term and an intraatomic exchange inter- 
action J term of multiband Hubbard-type minus a so- 
called double counting term, to subtract the electron- 
electron interaction already included in the LSDA . Our 
calculations are performed using the relativistic version 
of the "around-mean-field" LSDA + U method (AMF- 
LSDA+U) 0, as implemented in the full-potential lin- 
earized augmented-plane-wave (FP-LAPW) basis 0] in 
which the spin-orbit coupling (SOC) is taken into account 

m 

Before considering in detail the effect of Coulomb cor- 
relations for the Pu-Am alloys, we need to make sure 
that the AMF-LSDA+U is capable to deal with the 
end-points, namely <5-Pu and fee-Am. Applied to <5-Pu 
0, this approach accounted for all basic characteristics: 
equilibrium volume, bulk modulus, 5/-manifold binding 
energy, including the disinclination to magnetic order- 
ing, appearing due to a combined effect of SOC and 
Coulomb U, which lead to a nearly closed 5/ 5 / 2 sub- 
shell, hybridized with conduction electron states. We 
note that the non-magnetic 5-Pu can be obtained in fully- 
localized- limit (FLL)-LSDA+U calculations when no 
spin-polarization is allowed in the LSDA part and un- 
reasonably small value of the exchange J < 0.48 eV is 
used. 

Here we focus on the issue how can the same method 
describe similar (in terms of the 5/ occupancy), but pre- 
sumably much more localized 5/-states in Am metal. 
We restricted ourselves to the fee phase, which covers 
a broad concentration range from <5-Pu nearly to pure 
Am 0]. Also, it is the Am phase stable at elevated tem- 
peratures and at high pressure 0|. The Am atom has 
six / electrons with full j = 5/2 shell and zero magnetic 
moment. When the 5/ states are treated as non-spin- 
polarizcd valence states in GGA, the calculated equilib- 
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TABLE I: Ground state properties of /cc-Am. Given are 
the 5/ occupation number n^j , the spin moment M s , the 
orbital moment Mi, and the total magnetic moment Mj 
(spin+orbital ), and equilibrium volume (V eq ) and bulk mod- 
ulus {B). The experimental value of B, taken from Ref. l2l| . 
belongs to the dhcp structure. 
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rium volume V eq ~ 120 (a.u.) 3 j2£j is much lower (almost 
by 40 %) than the experimental value of 198 (a.u.) 3 |2lj . 
When spin polarization is included in GGA calculations 
|l2j . au gme nted by well known orbital polarization cor- 
rection^ (GGA+OP), the V eq = 170 (a.u.) 3 improves 
significantly. However, large spin M s = 6 fiB and orbital 
Mi = —1 /is moments are induced at the Am atom. This 
strong magnetic state, either ordered or disordered, was 
never seen experimentally. 

Another line to improve V eq assumes the 5/-manifold 
to be a part of an atomic-core (the "open-core" approx- 
imation), and neglects an overlap between the 5 /-states 
and the rest of valence electrons. It increases the vol- 
ume to 180 [13 or 202 [H (a.u.) 3 in agreement with 
the experiment. While the "open-core" approximation is 
based on a very strong constraint and corresponds to a 
non-physical limit of infinitely large Coulomb repulsion in 
the f-shell, it illustrates the impact of the 5/-localization 
for Am. Alternatively, self-interaction correction (SIC) 
to LSDA can be employed to localize the 5 .f 6 m anifold. 
The LSDA-SIC yields the V eq of 213 (a.u.) 3 jH which is 
somewhat larger than experimental value. Similarly to 
the "open-core" , the LSDA-SIC gives a localized Am 5/ 6 
shell. While being ab-initio and working well at Am end, 
the LSDA-SIC strongly overestimates the localization in 
S-Pu yielding V eq = 218 (a.u.) 3 [23| exceeding the experi- 
mental value by about 30 %. This illustrates the difficulty 
of LSDA-SIC method to deal with electron-electron cor- 
relations in <5-Pu, which limits further the applicability 
of LSDA-SIC to Pu-Am alloys. 

We performed AMF-LSDA+U calculations for Am, 
starting from the strongly spin-polarized LSDA charge 
and spin densities and on-site spin and orbital occu- 
pations, and without any constraints. We varied the 
Coulomb U from 3 eV to 4 eV, and chose the intra-atomic 
exchange parameter J=0.75 eV, i.e. in the range of com- 
monly accepted values for Am. The calculations con- 
verged to practically zero magnetic moment with residual 
M s and |M/| less than 0.01 \ib, and 5/-manifold occupa- 



tion n 5 f = 5.9. Also, V eq = 186 (a.u.) 3 and bulk modulus 
B = 55.1 GPa are calculated for U = 4 eV, in a reason- 
able agreement (considering the different structure type) 
with experiment (see Table I.). 
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FIG. 1: Color-on-line, (top) The DOS (TDOS-total and 
fDOS-projected) for /cc-Am from LSDA (inset), and AMF- 
LSDA+U= 4 eV calculations, (bottom) The TDOS and fDOS 
for PusAm alloy 



In Fig. 1 (top) we show the total (TDOS) and /-states 
projected densities of states (fDOS) resulting from the 
/cc-Am LSDA and AMF-LSDA+U=4 eV calculations. 
The strong spin-polarization found in LSDA disappears 
in AMF-LSDA+U, resulting in a truly non-magnetic 
(S z = L z = J z = 0) 5/ 6 ground state. The 5/ manifold 
binding energy is calculated as « 3.5 eV (for U = 3 eV), 
4 eV (U = 4 eV). The LSDA+U method is not based on 
any kind of atomic coupling scheme (LS or jj), it rather 
determines the set of single-particle orbitals that mini- 
mize variationally the total energy. The AMF-LSDA+U 
result can be interpreted as yielding the J = singlet 
ground-state configuration that corresponds to the jj- 
coupled Slater determinant formed of fully and equally 
populated six j = 5/2 orbitals. 

Although the DOS cannot be directly compared with 
5/ features in valcncc-band photoelectron spectra (PES), 
as - for localized 5/ states - those reflect not the 5/™ 
ground state, but the whole 5 f"" 1 multiplet, the energy 
range is undoubtedly correct pi. l2fij . Comparing to re- 
sults of Pu calculations (see Fig. 3 of Ref. p), the Am 
5/ states are not only at higher binding energy, but also 
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conspicuously narrower, with the width not exceeding 0.1 
eV. The results obtained for Am at ambient pressure can 
be definitely characterized as 5/-localized. Shift of the 
5/ states in AMF-LSDA+U down from the vicinity of 
the Fermi energy (Ep) decreases their contribution to 
the chemical bonding and increases V eq . 

The situation of the 5/ states changes dramatically 
when the volume is reduced from 198 (a.u.) 3 to 133 
(a.u.) 3 , which corresponds to the volume of the high pres- 
sure phase AmIV. This phase is generally assumed, as the 
orthorhombic structure reminiscent of a-U, to have itin- 
erant 5/ states Indeed, our calculations exhibit the 
5/5/2 band broadened nearly up to 1 eV. The width of the 
5/5/2 band becomes thus comparable to (5-Pu, although 
there remains a large difference in binding energies: in 
5-Pu, the 5/5/2 band is located around 1 eV below Ep, 
while in Am it is centered at 3.5 eV even for the high- 
pressure phase, just below the lower edge of the valence 
band. 

Next, we turn to the salient aspect of our investi- 
gation, the AMF-LSDA+U calculations of Pua,Ami_ x 
(x = 0.25, 0.5, 0.75) alloys. In all calculations, the 
Coulomb U= 4 eV, and exchange J=0.75 eV for Am 
atom, and £7=4 eV, J=0.70 eV for Pu atom were used. 
The /cc-based super-cells were used, LI2 for PusAm 
(Puo.75Amo.25) and PuAm 3 (Puo.25Amo.75), and Ll 
for PuAm (Puo.5Am.05). The lattice parameters, 8.905 
(a.u.) 3 for Puo.75Amo.25, 9.018 (a.u.) 3 for Puo.sAmn 
9.132 (a.u.) 3 for Puo.25Amo.75, were taken from Ref. [l| 

The LSDA calculations produce magnetic solutions for 
all considered Pu-Am alloys, similar to those reported 
in Ref. 0- Starting from these spin-polarized LSDA 
charge and spin densities and on-site spin and orbital 
occupations, we apply AMF-LSDA+U. The magnetism 
collapses in all cases, and Pu-Am alloys become non- 
magnetic, similarly to fcc-Am and <5-Pu. The calculated 
DOS (total per unit cell and 5/-projected (fDOS)) are 
shown in Fig.l (bottom) for Puo.75Amo.25 (PU3A111) al- 
loy. The Am 5 /-manifold is located at ~4 eV below 
E F , similarly to AMF-LSDA+U DOS for fee- Am (see 
Fig. l(top)). The Pu-5/ states are located in a re- 
gion from -2 eV to 0.5 eV and peaked at ~1 eV be- 
low Ep, similarly to the DOS of 5-Pu shown in Fig. 
3 of Ref. H (one also should keep in mind that there 
are small DOS changes due to the change in the lattice 
parameter). Importantly, the alloy DOS is almost ex- 
actly a weighted superposition of pure fee- Am and 5-Pu 
DOS showing that there is a very little interaction be- 
tween Am- and Pu-5/ manifolds. The calculated DOS 
for PU0.5A1110.5 and Puo.25Amo.75 look qualitatively very 
similar to Puo.75Amo.25 DOS shown in Fig. 1 (bottom), 
with well separated Am and Pu 5/-manifolds, whose rel- 
ative heights arc proportional to the Am and Pu atom 
concentration, again supporting the smallness of the in- 
teraction between the Am- and Pu-5/ states. The differ- 
ence between equilibrium volumes of S-Pu and fcc-Am 
is relatively small (« 15 %) and no significant deviation 
from Vcgard's law is found in experiment |l8j| . Consid- 
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FIG. 2: Color-on-line. Upper panel:valence-band photoelec- 
tron spectra (40.81 eV photon energy) for Am and <5-Pu. The 
lower panel displays the calculated ground-state DOS for Am 
and 5-Pu and the position of excitations obtained from the 
atomic calculations of Am (black). 



cring the lattice compression with a decrease of Am con- 
tent, we did not find any sizeable effect on Am-5/ states. 
The lattice compression is simply too small compared to 
the changes at hight pressure, and Am remains still prac- 
tically in the Ami phase, i.e. below the first phase transi- 
tion at 6 GPa. Similarly, there is no significant change in 
Pu-5/ states found with a decrease of Am content, just 
a small (« 0.2 eV) decrease in the binding energy. It is 
also noteworthy that we did not find any sizable change 
in the Pu-5/ and Am-5/ manifold occupations (715/) |26j |. 

In Fig. 2 we show the valence-band PES of 5-Pu, Am, 
and Puo.07Amo.33, obtained with the photoexcitation en- 
ergy 40.81 eV, emphasizing the emission from the 5/ 
states. Although in the case of band systems PES pro- 
vides practically a snapshot of the density of one-electron 
states, final-state effects are important for localized as 
well as narrow-band systems. In analogy to rare-earths, 
the valence-band spectrum of Am should exhibit the mul- 
tiplct of excited states. To address these effects, the one- 
particle Green's function and corresponding spectral den- 
sity were calculated for Am atom. In these calculations 
we used exact diagonalization of the many-body hamil- 
tonian for /"-shell closely following Ref. [28j , adding the 
SOC (£=0.346 eV), and using the values U= 4 eV and 
J= 0.75 eV, as in the AMF-LSDA+U calculations. Also, 
fcgT = 0.1 eV and (Nf) = 6 were used to determine the 
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chemical potential (Ep). The calculated spectral density 
shown in Fig. 2 (Am - photoexcitations) and a gree s very 
well with similar recent calculations of Svane [23. It is 
in a conspicuously good agreement with experimental PE 
spectra. Such comparison cannot be done in the case of 
Pu, as the use of the atomic limit for the Pu 5/ states 
is not justified. The Am-/ 6 many-body ground state is 
found to have J = character with both S = L = 2.471, 
while neither S nor L is a good quantum number. It dif- 
fers from the the result of Ref. |3(| that yields LS-couplcd 
7 F ground state. 

Although the experimental PES cannot be accounted 
for by the LSDA+U theory, one can see that the calcu- 
lated spectrum of the Pu-Am alloy corresponds to the 
weighted average of spectra of pure Pu and Am, sim- 
ilarly to the PE spectra. As shown recently J27J, the 
Pu-4/ spectra of Pu-Am alloys remain practically identi- 
cal irrespective of Am concentration at least up to « 35 
% Am. Also, very recent results of specific-heat studies 
Q indicate no enhancement of the 7 coefficient of the 
electronic specific heat. We can conclude that the lattice 
expansion caused by Am is not driving the Pu 5/ states 
to further localization and no magnetic moments on Pu 
in 5-Pu-Am alloys are induced thereby. Both from the 
experiment and theory, we deduce that the character of 
the Pu-5/ states does not change with the Am doping. 
The experimental PES features close to Ep can not be 
explained in terms of simple band states, similarly to the 
case of S-Pu for which the dynamical effects become very 
important pj. 



To summarize, we demonstrated that the around- 
mean-field version of the LSDA+U method gives a uni- 
fied picture of the electronic structure of Am, S-Pu and 
(5-Pu-Am alloys. For fcc-Am our calculations yield a 
non-magnetic ground state which is almost exactly a lo- 
calized f e J — singlet state, with equilibrium volume 
and bulk modulus in a good agreement with experimental 
data. Based on the calculations, we expect that this non- 
magnetic state will be preserved even for the pressure- 
delocalized state of Am, in accord with electrical resis- 
tivity 0] lacking any clear anomaly of magnetic origin. 

Neither magnetism nor further localization of the Pu- 
5/ states is found in fee Pu-Am alloys on the LSDA+U 
level, i.e. without introducing any additional effects 
like magnetic moment disorder or magnetic fluctuations, 
which are disproved recently by /xSR experiments |3lj . 
We are aware that our static mean-field approach does 
not account for a possible Kondo screening of the mag- 
netic moment fluctuations, which are important in other 
f-clcctron systems (e.g. a-Ce). 
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